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Abstract. Two-dimensional directed site percolation is studied in systems directed 
along the rr— axis and limited by a free surface at y—±Cx''. Scaling considerations 
show that the surface is a relevant perturbation to the local critical behaviour 
when ft<i/z where z=i^||/i^ is the dynamical exponent. The tip— to-bulk order 
parameter correlation function is calculated in the mean-field approximation. The 
tip percolation probability and the fractal dimensions of critical clusters are obtained 
through Monte-Carlo simulations. The tip order parameter has a nonuniversal, C— 
dependent, scaling dimension in the marginal case, k—l/z, and displays a stretched 
exponential behaviour when the perturbation is relevant. The fe-dependence of the 
fractal dimensions in the relevant case is in agreement with the results of a blob 
picture approach. 



1. Introduction 

The shape of the free surface limiting a system may influence its local critical behaviour 
at a bulk second order phase transition, provided the deviation from the flat surface 
is long range. For corners in two dimensions, wedges or cones in three dimensions, a 
marginal behaviour is obtained, with local critical exponents depending on the opening 
angle, in isotropic critical systems [1-9]. This result is linked to the invariance of the 
shapes under isotropic rescaling (see [10] for a review). 

With parabolic shapes, the surface introduces a relevant perturbation to the flat 
surface fixed point leading, locally, to stretched exponential behaviour for the order 
parameter and the correlation functions [11]. The Ising model [11-13], the self- 
avoiding-walk [14] and ordinary percolation [15] have been studied for this geometry 
in two dimensions. 

In anisotropic systems the correlation length diverges as ^n ^ t^'^n along a time- 
like direction and as ^j^^t^" in the transverse directions with a dynamical exponent 
z = Va jv [16] . Covariance under a change of the length scales then requires two different 
scaling factors, 611 —\f and h^ — h [17]. In this way, the relation between the correlation 
lengths fii ^ (Cx)^ is preserved. As a consequence, scale invariant shapes are different 
from the isotropic case: they now correspond to parabolic-like surfaces [18]. For 
example, the marginal shape for the directed walk, with z = 2, is the true parabola in 
two dimensions, a paraboloid or a parabolic wedge in three dimensions [18,19]. 
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In the present work, we study the two-dimensional directed site percolation 
problem inside a parabolic-like system. The scaling behaviour is discussed in section 2. 
The problem is solved in a mean-field approximation in section 3. The results of 
Monte-Carlo simulations are presented in section 4 and we end with a discussion in 
section 5. 



2. Scaling considerations 

We consider a system displaying anisotropic critical behaviour and limited by a 
free surface a.t y ~ ztCx'^ in the {x, y)-plane where x is the time-like coordinate. 
Under rescaling, with x' = x/b^ and y' = y/b, the geometrical constant C transforms 
according to 

^=&'-^'^- (2.1) 

When k > 1/z, 1/C renormalizes to zero, i. e. the system flows towards the flat surface 
geometry. The perturbation introduced by the free surface is irrelevant with respect to 
the flat surface fixed point. On the contrary, when k< 1/z, 1/C grows under rescaling 
and the system becomes locally narrower. The perturbation is then relevant and one 
expects a strong reduction of the order at the tip associated with a new type of local 
critical behaviour. The value k = 1/z corresponds to a scale-invariant shape leading 
to a marginal local critical behaviour with C-dependent exponents. For an isotropic 
system, z — l and the marginal shape, with fc=l, corresponds to the corner geometry 
mentioned above. 

Let ttiq be the tip order parameter, with scaling dimension a;^ at the flat surface 
fixed point, on a finite system with size L along the x-axis. Under a change of scale, 
it transforms as 

4- i) -"--»«('■'"'• ^'f) <") 

where t is the thermal scaling field and 1/C, which is vanishing at the reference fixed 
point like 1/L in finite-size scaling, is treated as a new scaling field introduced by the 
parabolic free surface. One may also notice that, for a directed problem, the system 
does not see a flat free surface at x = and x^^^ is also the bulk scaling dimension of 
the order parameter. 

Taking b = t^" ^ one obtains 

t,-,-\^tPmJ^,—\, l^^C^ni-zk) (2.3) 

'C l) ^ \t-''n ' t-^n J ' ^ ^ ' 

where (3 — vx^ and the shape of the system introduces a new length scale given by 
Iq as long as the shape is not scale-invariant, i. e. k^ 1/z. In the same way, taking 
b = lfj in equation (2.2), one obtains the following finite-size scaling behaviour at the 
critical point t = 0: 

«.h-l)-'c"'""'c(r)^ (2.4) 
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Figure 1. Directed site percolation on the square lattice: (a) The system is directed 
along the diagonal of the square lattice, (b-d) The site at (x + l,y) (full circle) is 
connected to the origin if at least one of its nearest neighbours at x is connected to 
the origin (full squares). 

In the marginal case, z = 1/fc, the scahng dimension of the tip order parameter 
becomes C-dependent and equation (2.2) is changed into 



m,[t,^,^)^b-^-(^^m^^^Jb^/'^t^' 



(2.5) 



The order parameter correlation function between the origin and a point at (x, y) 
transforms as 



G 



(,.4) ^. -«(.-. J, f,^) 



when L is infinite. With b — x'^, equation (2.6) leads to 



G{t,x,y,— ] =x 



^-2x^/z 



X y X 



t-^W x' Ir 



(2.6) 



(2.7) 



These scaling forms will be used later to analyze the mean-field results and the 
Monte-Carlo data. 



3. Mean field approximation 

We consider the site percolation problem on a square lattice, directed along a 
diagonal of the squares as shown in figure la, with a site occupation probability p. 
The order parameter correlation function is the probability density P(x, y) for a site 
at (x, y) to be connected to the origin. The site at (a;+l, y) will be connected too if it is 
occupied and if at least one of its two nearest neighbours at x are themselves connected 
to the origin (figures Ib-d). In a mean-field approximation [20] where correlations 
between the Ps on different sites are neglected, this occurs with probability 



P 



(3.1) 



{x + l,y) ^p[p{x, y + l)[l- P{x, y - 1)] + P{x, y - l)[l - P{x, y + 1)] 
+ P{x,y + l)P{x,y-l)^ 
where the different terms on the right correspond to the last three diagrams in figure 1. 
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Going to the continuum limit and neglecting nonlinear terms involving derivatives, 
one is led to the following differential equation: 

8P 8^P 

^-P^ + (2P-1)P-P^^ (3.2) 

The order parameter is the probability Pq for the site at the origin to belong to an 
infinite cluster. For the un-confined system, it is given by the homogeneous stationary 
solution of equation (3.2) as 

P„= ^ ^-^^ 2 (3.3) 

L p<Pc 

and vanishes at the percolation threshold p^ with an exponent /3'"-' =f. 

Below the threshold and far from the origin, P{x, y) <S^ 1 and equation (3.2) can 
be linearized, leading to: 

dP 1 d^P 

-^-— + 2(p-pJP^ (3.4) 

Through the change of function P{x, y) -^ exp [2(P|,— p)x] P{x, y) an ordinary diffusion 
equation is obtained, so that 

Pix, y) = e-2(Pc-p)- """^A^) (35) 

for the un confined system. In d+l dimensions the power-law decay would be changed 
into x^'^l'^ . From a comparison with equation (2.7) where t is now |p— p^,], one deduces 
the following exponents for directed percolation, 

z/^^ = 1 v^^ = 1/2 z^^ ^2 xZ^ ^'^ (3.6) 

where v"^^ follows from scaling. The scaling law (3 = vx^^ is verified for the 
mean-field (Gaussian) exponents only at the upper critical dimension d^ + 1 which, 
according to (3.3) and (3.6), is equal to 5 for directed percolation. 

On a parabolic system, we use the new variables x and (^{x, y) = y/x^ for which 
the free surface is shifted to (l = ±C and equation (3.4) is changed into 

dP 1 d'^P CdP , , „ , , 

with the boundary condition P{x, Q = ±C) = 0. Through the change of function 

P{x, = exp [ - 2(p, - p)x - ^ex:"^-'] Q{x, Q) (3.8) 

equation (3.7) leads to 
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for which the variables separate when k = 1, 1/2 or 0. According to equation (2.1), 
in mean-field, these values of k just correspond to irrelevant, marginal and relevant 
perturbations. 

For k = l, i. e. in the corner geometry, (3.9) gives 



,dQ X 



Id^Q 



X 



p — Q-2(Pc-p)2:-J/^/2a;Q 



dx 2^ 2 dC 
and the correlation function, which is even in y, takes the form 



(3.10) 



P(x, y) = e-2(p<=-p)=^ 



^2'ITX 



E^« 



exp 



(2n + l)27r2 



8C2x 



(2n + l)7ry 



2Cx 



(3.11) 



This is just the un-confined solution (3.5), modulated by a function which depends 
on the two last variables v and w of the scaling function g{u,v,w) in equation (2.7), 
with here z = 2 and l^^C^^. The critical behaviour is the same as for un-confined 
percolation as expected for an irrelevant perturbation. 

For the true parabola which is the marginal geometry, one may use equation (3.7) 
with fc = l/2 to obtain 



dQld^ (dQ 



c 



y 

pl/2 



P 



-2{pa-p)x 



Q 



(3.12) 



which is of the form studied in [18] for the directed walk problem. Writing Q{x, () - 
(j){x)i}'{(^) leads to the following eigenvalue problem for ipiC) 



2de+2dC^"^^ 



(3.13) 



with (f){x)^x ^ . The solution, which is even and regular at the origin, can be written 
as a series expansion, X^^o'^nC^"- According to (3.13), the coefficients satisfy 



(A2 + n) 



'"n+l 



(n+l)(2n+l) 



(3.14) 



which is the recursion relation for the coefficients of the degenerate hypergeometric 
function 



iJ^i 



1 ^2 



^^ ' 1.3. •• -(271-1) 



) C 



2n 



(3.15) 



The boundary condition '0(C') — gives the eigenvalues A„ which are the zeros of 
1-^1 [-^n' 1/2; —^2/2] . The solution is obtained as the eigenvalue expansion 



P(a;,y)=e-2(P-P)-^B„ 



1^1 



n=0 



X2 



2x 



(3.16) 



The behaviour at large x is governed by the first term in this expansion which decays 
as x~'^o, i. e. with a C-dependent exponent as expected for a marginal perturbation. 
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Figure 2. Scaling dimension of the tip order parameter as a function of C in the 
mean-field approximation when the shape is marginal. The local exponent xJJJ^ (C) 
diverges when C— >0 and goes to its unperturbed value x^^ = 1/2 when C— >oo. 



The dimension of the tip-to-bulk correlation function is the sum of the tip and bulk 
order parameter dimensions, the first one being variable. Comparing with the form 



of the decay in (2.7) gives Aq = 
parameter dimension is given by 



mf 



r-mf 



{0)^2X1-1 



{C)+x^^]/z"^f and, using (3.6), the tip order 



(3.17) 



Its dependence on C is shown in figure 2. 

Analytical results can be obtained only in limiting cases which have already been 
discussed in [18]. When C is infinite, Aq = 1/2, only the first term in the expansion 
remains, which satisfies the initial and boundary conditions, giving back the free 
solution in equation (3.5) since iF-^{^, i; — j/^/2x) = exp(— j/^/2x). For large C-values, 
the tip exponent takes the following form 



„mf(^x _ 



1 



(C) = T. + 



l^M-^ 



[l + 0(e)] 



(3.18) 



where s is the correction term itself. For narrow systems, the hypcrgcomctric function 
gives a cosine to leading order in C^. One obtains the following asymptotic behaviour 
in X 



P{x,y)^x-^"/'^'" 



cos 



\2CV^J 



(3.19) 



and the tip exponent diverges as 7r^/4C^. 

For the strip geometry, which corresponds to a relevant perturbation, one can 
use (3.9) with fc = or, more simply, go back to the original equation (3.4) giving 



dQ 
dx 



2 92/2 



p 



'2{pa-p)x 



Q 



(3.20) 
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The solution satisfying the initial condition P{0,y) — 6{y) reads 

Pix,y)^ C ^Z^^^P[ ^2 Jc°s[ ^ J (3.21) 

n=0 

in agreement with the scaling form (2.7) for the mean-field exponents. The 
perturbation due to the free surface induces a simple exponential decay of the 
correlations at p^. Actually, the system is now one-dimensional and this term 
corresponds to a shift of the percolation threshold by (tt/AC)'^. 

For < fc < 1/2, the simple exponential is expected to be replaced by a stretched 
one, involving some power of the last scaling variable w = x/l(j in equation (2.7). The 
power goes to 1 when fc = and has to vanish when the system is marginal (A: = 1/2) 
in order to give a C-dependent decay exponent. A good candidate for the asymptotic 
behaviour is [18] 

D/ ^ exp[-2(p^-p)a;] / tt^ a;i-2fe \ f Try \ 

It has the proper scaling form (2.7) with mean-field exponents and interpolates 
between the strip result (3.21) and the small-C result (3.19) for the parabola. 



4. Monte Carlo simulations 

As in the mean-field approach of the last section, the Monte-Carlo (MC) 
simulations were performed for site percolation on a square lattice along the diagonal 
direction. Accurate estimates of the percolation threshold and directed percolation 
exponents have been deduced from series expansion in [21] with 

p^ = 0.705489(4) v ^ 1.097(2) vn = 1.734(2) 

(4.1) 
7 = 2.278(2) 2 = 1.581(3) /3 = 0.276(3) 

where the last two values follow from scaling laws. 

On the lattice, the sites inside or on the curve y = ±Cx'' are considered to belong 
to the system. The tip order parameter was calculated on parabolic-like systems with 
size L along the a;-axis. On these finite-size systems, Pq is the probability that a 
cluster grown from an occupied site at x = y = reaches x~L. Such a cluster is called 
"infinite". "Finite" clusters are those dying before x = L. 

In order to spare computer time the MC-algorithm ensures that a cluster reaches 
a; = L as fast as possible. This is done by always trying to occupy those sites with the 
largest x-coordinate. The coordinates of sites from which the growth may eventually 
continue but which are closer to the origin are stored into a list of active sites. The 
growth stops once the size of the cluster is equal to L, even if it still contains active 
sites. Thus, the complete infinite cluster has not been built. Subsequently, if there are 
no active sites left, a complete finite cluster has been grown. An example is shown in 
figure 3. 
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Figure 3. Directed percolation cluster generated using the Monte-Carlo algorithm 
described in the text. Here, L = 1000, fc = l/2 and C = 4. 



0.8 - 



0.6 



0.4 



0.2 - 



— 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 






Pc 


— 1 — m — 1 — 


1 1 1 


-T — 1 — 1 — \ — r- 


-1 — I— 

Rill- 


— 










.»'- 


f, ^„»° 


— 


C = 4 










t'''**'. 


"'' 


^pD- 


C = 3.5 








«:•: 


* ^ » " 


^o" 


° _ 


C=3 










^.•° 


□ °° 




C = 2.5 








s"^ 4 




□ ° 


- 


C = 2 








'•*** ."* 


^„ 




- 


- C=1.S 

C = l \ 


7 




:■: 


' .' 


_."°° 




D^"" 


- 








" 




^„°"° 


- 


- k = l/2 




r- 


• 


□ ° 


^.o"" 




" 




h 


^ 




= ?°° 1 


1 1 







0.65 



0.7 



0.75 



0.8 



Figure 4. Tip percolation probability on a parabola (fc = l/2) with length L = 1000 
as a function of p. The Monte— Carlo results are averaged over 10^ samples for 7 values 
of C between 4 and 1 from top to bottom. The non-vanishing percolation probability 
at the percolation threshold (dashed line) indicates strong finite-size effects which 
are increasing with C. 



With this algorithm, which is a modified version of the one used in [15] for ordinary 
percolation, the computation time for increasing p > p^ remains almost constant. Thus 
we could compute Pq as a function of p up to p = 0.8 (A; = 1/2, 1/z) for L = 1000 taking 
averages over 10''' samples. The results are shown in figures 4 and 5 for values of C 
ranging from 1 to 4. 



Due to the anisotropy, finite-size effects remain important. This appears in the 
non-vanishing values of Pq{Pc) and increases with C, i. e. when the system opens. 
According to the scaling relation (2.3), P^t"^ is an universal function of lc/^\\ and 
i/^l|. This has been verified in the relevant case k = 1/2 for the first of these 
scaling variables. A good data collapse is obtained in figure 6 for C — 1, 1.5, ... ,4 
and L = 50, 100, . . . , 1000 with a fixed ratio L/^n =0.34. 

In order to check the expected stretched exponential behaviour in the relevant 
case, we performed a finite-size scaling study at the critical point. The tip percolation 
probability, for A; =1/2 and 1<C<4, was calculated at the percolation threshold on 
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Figure 5. As in figure 4 for the marginal case, k = l/z. 
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Figure 6. Scaling function Pot ^ versus Ic /ill in the relevant case (fc = 1/2). A good 
data collapse is obtained for the different values of C (same symbols as in figure 4). 



systems with size 50<L< 1000, taking averages over 10^ MC samples. According to 
equation (2.4), P^V^ is a function of L/l^^ which depends only on k. The scaling 
behaviour is shown in figure 7 where the deviations at small C-values is likely to be due 
to corrections to scaling. Actually, the reference fixed point for which equation (2.4) 
was written corresponds to C— >oo. In this semi-logarithmic plot, a linear dependence 
on {L/l(jY~^^ is obtained. It corresponds to a stretched exponential behaviour for the 
scaling function m(j{L/l(j) although some power in front of the exponential cannot be 
excluded. Other simple combinations of z and k for the stretching exponent did not 
lead to a linear behaviour. 

In the marginal case, k — 1/z, the scaling dimension of the tip order parameter 
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Figure 7. Finite-size scaling of the tip order parameter at the eritieal point for a 
relevant surface shape (A; = 1/2, same symbols as in figure 4 for C). The almost linear 
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, function of (L/Iq)^ '^^ indicates a stretched exponential 



behaviour of the tip percolation probability. 
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Figure 8. Scaling dimension of the tip order parameter as a function of C obtained 
through finite— size scaling when the shape is marginal {k = l/z). The local exponent 
Xm(C) diverges when C^O and goes to its unperturbed value Xm =0.252(3) when 



is expected to vary with C. Equation (2.5), with 6= L^/^, leads to the finite-size 
behaviour Pq(L)~L~^"'('-^-''^ at p = Pc which can be used to determine x^{C). This 
has been done taking averages for Pg (L) over 2.10^ samples with L — 5Q0, 520, . . . , 1000. 
The tip exponent was deduced from the asymptotic slopes of log-log plots. The results 
are shown in figure 8. The dependence on C is qualitatively similar to the mean-field 
one in figure 2: the scaling dimension diverges when C vanishes and it decreases to 
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Figure 9. Fractal dimensions versus k for finite critical percolation clusters starting 
at the tip. The variation in the relevant regime, fc< 1/2: = 0.633(2), is vifcU described 
by the blob picture approach (lines) when C decreases. The deviations for large 
C-values are due to the finite size of the system: only a few blobs can develop and 
the fractal dimensions are then closer to the unperturbed ones. 



the bulk value x^—(3/iy — 0.252{3) when C goes to infinity. 

The fractal structure of finite critical clusters was also studied. Since details can 
be found in reference [22], here we only give a brief summary of our results. As 
in figure 3, 2.10^ percolation clusters starting from the tip were generated at the 
percolation threshold, inside a parabolic system with size L=1000 for C = 2,3,4 and 
values of k ranging from 0.25 to 0.75. The mean square radii of gyration, {X'^) in the 



in the transverse direction, were determined for s-site "finite" 



x-direction and (Y^) 
clusters, i. e. 

For an unconfined system, the cluster size behaves as 



clusters with x^^^ < L. 



X^ 



K 



(4.2) 



where X^ = (X^)^/^, Y^ = (Y^y^^. The exponents are fractal dimensions which, 
extending an argument of Stauffer [23], are related to the directed percolation 
exponents through 



11 



/3 + 7 
^11 



1.473(2) 



/3 + 7 



2.329(3). 



(4.3) 



The second value is greater than the dimension of the system due to the anisotropy as 
explained in reference [22]. A single fractal dimension df may be also defined using the 
characteristic length I ^ ('||^j_)^^^ ^ s^^'^f associated with the surface of the cluster 
so that [20] 2/df^l/d^^ + 1/d^ and then, d^ = 1.805(2) < 2. 

Inside a parabolic system, one expects a similar behaviour with fc-dependent 
fractal dimensions dJk), dj^{k), when the shape is a relevant perturbation. Log-log 
plots of the MC results are linear for cluster sizes between 2^ and 2^^-2^'^. Deviations 
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at larger sizes are due to finite-size effects. Tire fractal dimensions are shown in 
figure 9 as functions of k. The form of the variation in the relevant regime, k < 1/z, 
can be obtained using a blob picture approach [24-26,14]. The cluster configuration 
is supposed to follow from the piling up of anisotropic blobs inside the parabolic 
system. Within each blob the structure is the same as for unconfined clusters with 
the unperturbed values of the fractal dimensions. In this way one obtains [22] 

dii(fc) = l + zfc(d|| -1), rf|(fc) = ^^, 0<k<l/z (4.4) 

II II fc 

in good agreement with the Monte-Carlo data as shown in figure 9. 



5. Conclusion 

The local critical behaviour for directed percolation has been investigated at the 
tip of two-dimensional parabolic-like systems. The problem has been treated using 
mean-field theory and MC-simulations. In the marginal case, k= 1/ z, the percolation 
probability and the order parameter correlation function display a non-universal local 
critical behaviour: the local exponents vary continuously with the shape of the system. 
In the relevant case, k < 1/z, the critical tip percolation probability is a stretched 
exponential function of the reduced size L/l^j and finite critical clusters have k— 
dependent fractal dimensions. 

Although we were unable, due to strong finite-size effects, to check a similar 
stretched exponential dependence on £.i\/lc in the relevant case, the following general 
picture seems to emerge from our mean-field and MC results, exact and MC results 
for the Ising model [10-13] and ordinary percolation [15]: 

• In isotropic critical systems, conformal transformations, from the half-space 
to the parabola, of the correlation function and the order parameter profile 
(with appropriate fixed boundary conditions) show that these quantities become 
stretched exponential functions of x/l(j or L/l^-, with a stretching exponent 1 — /c. 
According to our mean-field and MC-results, this exponent should be changed 
into 1 — zk in anisotropic critical systems. Conformal methods can no longer be 
used with anisotropic critical systems for which the conformal group is replaced 
by the Schrodinger group (when z = 2) [27] . 

• In isotropic off-critical systems, the tip order parameter is a stretched exponential 
function oi ^/Iq (with there l^j — C^^'^^^''^) with a stretching exponent {l~k)/k. 
This follows from exact results on the Ising model [11,12] and also from an heuristic 
argument assuming that the tip order, which is induced by the bulk order at a 
distance D ^ {i/CY^^ (where the width of the system is of the order of the 
bulk correlation length), decays with D like the critical tip-to-bulk correlation 
function [10,11]. The same argument, applied to the anisotropic off-critical 
system, gives a stretched exponential function of ^\\/lc with a stretching exponent 
{\ — zk)/zk. Such a behaviour could not be checked numerically here, since our 
finite off-critical system involves some unknown combination of the two scaling 
variables L/l(j and Cii/'c- 

Further numerical and analytical work on anisotropic systems would be useful to 
confirm the last point. 
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